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We theoretically study the spin-polarized transport through a single-molecule magnet, which 
is weakly coupled to ferromagnetic leads, by means of the rate-equation approach. We consider 
both the ferromagnetic and antiferromagnetic exchange-couplings between the molecular magnet 
and transported electron-spin in the nonlinear tunneling regime. For the ferromagnetic exchange- 
coupling, spin current exhibits step- and basin-like behaviors in the parallel and antiparallel con- 
figurations respectively. An interesting observation is that the polarization reversal of spin-current 
can be realized and manipulated by the variation of bias voltage in the case of antiferromagnetic 
exchange-coupling with antiparallel lead-configuration, which may be useful in the development of 
spintronic devices, while the bias voltage can only affect the magnitude of spin-polarization in the 
ferromagnetic coupling. 

PACS numbers: 75.50.Xx, 73.23.-b, 72.25.-b, 85.75.ci 



I. INTRODUCTION 

In the past few years, electron transport through mag- 
netic molecules, especially the single-molecule magnet 
(SMM), was intensively studied in both experimental 0- 
5] and theoretical [6H23I] aspects, which is stimulated by 
the fundamental importance as well as potential ap- 
plications in molecular spintronics. Many fascinating 
properties have been found, such as complex tunneling 
spectra [(| 0], negative differential conductance (NDC)[l[ 
l8l-[lfj{. Kondo effec t I13 . fl3j . Berry phase blockad e 1141 . full 
counting statistics [l5|, colossal spin fluctuations [l6j. and 
so on, which are resulted from the large spin-number and 
high anisotropy of SMM. The complete current suppres- 
sion and negative differential conductance are also con- 
firmed by the single-molecule-transistor measurement, 
which as a function of bias, gate voltage, and exter- 
nal magnetic field provide evidences of magnetic signa- 
tures of the SMM PI . Manipulation of spins in mag- 
netic molecules, which is based on the development of 
spin-controlling techniques, may result in new strategies 
of quantum-state control. It is demonstrated that the 
charge and spin states of molecule can be identified from 
the measured tunneling spectra due to the existence of 
exchange interaction between the local spin of magnetic 
molecule and the spin of tunneling electron, particularly 
the antiferromagnetic exchange interaction in a magnetic 
single- molecule transistor based on N@C60[3, @. On the 
other hand the tunneling current between magnetic elec- 
trodes can also control the orientation of molecular mag- 
net shown by the study of colossal spin fluctuations [la] . 

Based on the study of spin polarized transport through 
a SMM some new spintronic devices are proposed, which 
display many interesting effects. In particular the spin- 
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FIG. 1: (Color online) The schematic diagram of quantum 
transport through a SMM. 



diode behavior [Hj and tunnel magnetoresistance (TMR) 
[HI, [13] are shown to be resulted from the exchange cou- 
pling between the lowest-unoccupied-malecular-orbital 
(LUMO) level and the core spin. Furthermore a large 
negative TMR is predicted in the case of antiferromag- 
netic exchange- coupling [171 ]. The giant spin amplifica- 
tion and spin-blockade behavior are useful in molecu- 
lar spintronics to serve as a read-out mechanism [8|, |9| . 
It is also found that the magnetization of SMM can 
be reversed by spin polarized current jj| Il8l - [2fjj and 
spin-bias driven magnetization reversal is also observed 
theoretically [2l|. The most recent study is focused on the 
manipulation of SMM by thermal spin-transfer torque 
22]. Moreover a highly spin-polarized current can be 
generated by thermoelectric effects [23 1. 

To date, most researches are concentrated on the ma- 
nipulation of a SMM by the transported electron, which 
on the other hand may be used as a probe to explore 
the level structure of the SMM. The spin-polarized cur- 
rent through a SMM itself has received little atten- 
tion, while in the quantum-dot (QD) case it is exten- 
sively studied [24l-|26{ in terms of nonequilibrium Green's 
functions [25j and master-equation [2|| respectively. The 
QD acts as a spin-current diode giving rise to the spin 
blockade. The spin current in QD system shows quite 
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different characteristics compared with the charge trans- 
port [n, m Hi . 

We in the present paper study the spin-polarized trans- 
port through a SMM, which has more complex level- 
structure than the single QD. Moreover the exchange 
coupling between SMM and the transported electron-spin 
leads to additional dynamic mechanism to manipulate 
the spin-polarization of current. Both parallel and an- 
tiparallel configurations of lead magnetization are consid- 
ered along with the ferromagnetic and antiferromagnetic 
couplings. For the ferromagnetic exchange-coupling, the 
spin current variation with respect to the bias voltage 
exhibits a step-like curve in parallel configurations and 
basin-like behavior in the antiparallel configurations, re- 
spectively. The interesting observation in antiferromag- 
netic coupling case is that, the spin-polarization can be 
reversed with increasing bias voltage in antiparallel lead- 
configurations. 



II. EXCHANGE COUPLING AND 
TRANSITION RATE OF TUNNELING 

We consider a system which consists of a SMM coupled 
to two ferromagnetic metallic-electrodes (see Fig. 1.), 
which can be described by the Hamiltonian [H, lij, Il9j 



H — Hsmm + H, 



leads 



H 7 



(1) 



The first term in Eq. (1) concerning the SMM of easy- 
axis anisotropy with parameter K z > has the form 

Hsmm = ^{e-eV g )d) a d (T +Ud\d^d\d i - Js-S-K z {S z ) 2 , 

a 

(2) 

in which Js ■ S is the exchange interaction between elec- 
tron spin and the giant spin S of SMM with J being the 
exchange coupling parameter and s= J2 aa ' d\{a aa i /2)d a ' 
is the corresponding spin operator of electron (cr is the 
vector of Pauli matrices). d\{d a ) denotes the relevant 
electron creation (annihilation) operator and e is the 
single-electron energy of the LUMO level, which is tun- 
able by the gate voltage V g . U represents the Coulomb 
interaction of two electrons of opposite spins. The ex- 
change interaction can be of either ferromagnetic ( J > 0) 
or antiferromagnetic ( J < 0) type. 

The many-body states of electron-spin and molecule 
are expressed in terms of the eigenstates of operator 
Sj ot , n, Stot] m), where Stot—s + S denotes the total spin 
operator of electron and molecule, with n denoting the 
charge state of the SMM, Stot the total spin quantum- 
number and m the eigenvalues of Sf ot [J [H, [l?], fl9j - 
For the case of n — 0, the state |Q, Stot = S; m) = 
\®)orb ® \ m ) mo i an d the corresponding eigenenergy is 



'\0,S tot =S;m) 



= —Kzin 2 , while when 



2, the state is 



\2,S to t = S;m) = |tl) orh ® \™>) mol and the eigenenergy 
is s\ 2 ,s tot =S;m) = 2(e - eV g ) + U - K z m 2 . For n = 1 
and K z — J/2 > 0, the state and eigenenergy are found 
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FIG. 2: (Color online) Energy spectra of the SMM for ferro- 
magnetic (a) and antiferromagnetic (b) exchange- interactions 
as a function of magnetic quantum-number m with S = 2, 
K x = 0.05meV, U = lmeV, Vg = -0.3mV. 



It) 0r6 ® \m ~ V 2 )moi and £\l,S tot =S±l/2;m) = 



eV g ) + J/4 - K z (m 2 + 1/4) ± AE(m) respectively. For 
n = l and K z - J/2 < 0, we have \l,S to t = S =p 1/2; m) = 
42 \Dorb ® \m + l/2) mol + &g |t) or6 ® |m-l/2) moJ , 
and £|i,s tot =STi/2;m) = (e - eVg) + J/4 - Jf^rn 2 + 
1/4) ± AE(m), where A£(m) = [K Z (K Z - J)m 2 + 
(J/4) 2 (25+ 1) 2 ] 1/2 , and a^, &g> are effective Clebsch- 
Gordan coefficients [Tij. 

The SMM is weakly coupled to two ferromagnetic 
metallic leads with the Hamiltonian given by 



Hi 



I'. mis ^ ^ ^ ' ^ak^Qitj^akff; 

a— L,R ker 



where c^ k(T (c Q k<x) is the creation (annihilation) operator 
for an electron of spin-index a and wave vector k in the 
a lead. The spin polarization of ferromagnetic lead a is 
defined as P a = (p a+ - p a -)/(p a + + Pa-), with 
denoting the density of states for the majority (minority) 
electrons in the a lead. In this paper, we assume that 
the magnetization of ferromagnetic leads is collinear with 
the easy axis of SMM. The tunneling processes between 
the molecule and leads can be described by Hamiltonian 

H T = ^2[t a cl^d a + t* a dlc aka ], 

aker 

with t a denoting the tunnel matrix clement between 
the molecule and a lead. The spin-dependent tunnel 
coupling-strength is given by 

faff — 2np a a \ta\ : 

and the total coupling-strength is T a — J2 a T ao -. ^ n or " 
der to obtain transport properties in both the sequential 
and cotunneling regimes, we employ the rate-equation 
approach with the help of T-matrix [29l432| , which satis- 
fies the iterative equation 



T = Hi 



H 7 



Et-Hq 



-T. 



(3) 



as \l,S tot = S±l/2;m) 



i {±) 



\l) orb ® \m + l/2) .+ 



where Ej is the energy of initial state \I), and + denotes 
a small quantity in the retarded Green function. The 
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transition strength from initial state \I) to final state \F) 
can be evaluated by the perturbation expansion of T- 
matrix in terms of the generalized Fermi's golden rule[29j], 



T FI = ^\{F\T\I)\ 2 S(E F -E I ) 




_ 2tt 
~~ Ji 
5{E F - 



(F\H T 



Hi 



1 



Ei - H + i0+ 



-H 7 



(4) 



where |7) and \F) are actually product states of the 
electron-lead and molecule. After eliminating the lead 
degree of freedom, the sequential transition rate up to 
the second-order of tunneling Hamiltonian Hx can be 
obtained from the transition-strength formula 



H a )\{i'\rfM 



+ [l-f( Si - Si ,- tla )]\(i'\d a \i)\ 2 }, (5) 

which describes the transition of molecule from state \i) 
to \i'), due to spin-tr electron tunneling into or out of lead 
a. Where /i£ = — eV/2 , fin = eV/2, Si is the energy 
of state \i) and f(x) is the Fermi distribution function. 
For the SMM model we can find the sequential transition 
selection-rule : |An| = 1, \AS tot \ = 1/2, and \Am\ = 1/2. 
The fourth-order cotunneling transition rate 
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FIG. 3: (Color online) (a) Charge current I c , (b) differ- 
ential conductance G, (c) spin current I s , and (d) current 
polarization \ as a function of the bias voltage V for dif- 
ferent lead-polarizations p and ferromagnetic exchange cou- 
pling (J = 0.2meV) in case of parallel configuration of lead- 
magnetization. 



from the condition, ^f- = 0, with the transition rate 
together result in the current of spin- a through lead a: 



stands for the virtual transitions from molecular states 



+ (W l,i> , , - W*'? , )PA 

' \ vv acr.a'a' y ' a 1 a 1 ^aa) 1 l J ' 



(8) 



to while changing a spin-a electron of lead a into from which wg obtain ^ charge current ^ = ^ 



spin-cr' electron of lead a' . The divergence of above ex- 
pressions existed whenever energy denominators vanish 
has to be eliminated by some regularization procedures 
[3ll - l3ll | and in this paper we adopt the method described 
in Ref ,[36|. The cotunneling transition selection rules 
are seen to be |Ara| = 0, and |Am| = 0, ±1. 

For the weak coupling between lead and molecule, i.e. 
T aa <C ksT, electron transport can be considered as a 
stochastic Markovian process and thus the time evolution 
is described by the following rate equations 



at 



m 

dt 



= E -(w^^+w^)Pi 

act' ' aa' ' i' '^i 

+ E (W^ + W2,, a< r)Pi>, (7) 
aa' ' crcr' i' '^i 

where Pi is the probability of finding molecule in the 
many-body state i. The stationary probabilities obtained 



I a ± and the spin current I as = I a -\ — 7 Q j,. The current 
polarization is defined by \ — (I a -\ — lai )/(Iarf + lal)- 



ILL CURRENTS AND POLARIZATIONS IN 
THE NONLINEAR REGIME 

The numerical results of the charge current, differential 
conductance, spin current, and current polarization in 
the nonlinear regime are carried out with parameters cho- 
sen as S = 2, s d = O.bmeV, Vg = -0.3mV, J = 0.2meV, 
U = ImeV, K z = O.ObmeV and k B T = O.OimeV. We 
set the tunnel coupling between SMM and ferromagnetic 
leads to be T = T F = = O.OOlmey and assume 
the same polarizations of two leads i.e., pl = Pr = P- 
In addition, the current and differential conductance are 
scaled in units of Iq = 2eT/h and Go = 10~ 3 e 2 /h, re- 
spectively. With the chosen parameters the numerical 
energy-spectrum of SMM is obtained as in Fig. 2, from 
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FIG. 4: (Color online) (a) Charge current I c , (b) differen- 
tial conductance G, (c) spin current I s , and (d) current po- 
larization x as a function of the bias voltage V for differ- 
ent lead-polarizations p and antiferromagnetic exchange cou- 
pling (J = —0.2meV) in case of parallel configuration of lead- 
magnetization. 

which one can find that the energy-eigenvalue of SMM- 
states corresponding to the total spin- number Stot = 5/2 
is lower than that of Stot = 3/2 for the ferromagnetic 
exchange-interaction (Fig. 2(a)), while the situation is 
just opposite for the antiferromagnetic case (Fig. 2(b)). 
In what follows the spin currents are analyzed for two 
types of exchange interactions with ferromagnetic leads 
of both parallel and antiparallel configurations respec- 
tively. 



A. Parallel configuration 

For the ferromagnetic exchange-interaction, the bias- 
voltage dependences of charge current I c and differen- 
tial conductance G for different lead-polarizations p are 
shown in Fig. 3(a) and Fig. 3(b) respectively, where 
it is seen that the deviations of I c and G are negligi- 
bly small when the lead polarization varies in the re- 
gion p < 1. This observation is, as a matter of fact, the 
same as in the QD '27], where the resonant probability 
is independent of the lead-polarization p [37] . There ex- 
ist five evident resonant-peaks in Fig. 3(b), in which 
the first (corresponding to the transition |0,2;±2) <^> 
|1, 5/2; ±5/2)) and fifth ( |1, 5/2; ±5/2) & |2,2;±2)) 
peaks are higher than the other three resulted from the 
transitions: |0,2;±1) \1, 3/2; ±3/2) (2) , |0,2;±2) <S=> 



|1, 3/2; ±3/2) (3) and |1, 3/2; ±3/2) <S> |2,2;±2) (4), 
respectively [17j. However, for the full polarization p = 1, 
the current and differential conductance are quite dif- 
ferent since only one spin-up channel is involved in the 
transition without competition between majority and 
minority spin channels, which leads to dynamic spin- 
blockade .9, 38] . The bias- voltage dependence of the spin 
current I s and current-polarization \ f° r various lead- 
polarizations p are shown in Fig. 3(c) and Fig. 3(d). 
The spin current exhibits step-like behavior in parallel 
lead-configurations, which is also the same as in a single 
QD (27[, while with more steps due to the complex en- 
ergy levels of SMM. In addition, when lead-polarization 
p increases, the spin current varies more evidently along 
with an increasing magnitude of current polarization x 
in contrast to the charge current. The current polariza- 
tion x approaches the largest value at a very low bias 
voltage, which is dominated by elastic co-tunneling pro- 
cesses through majority-majority and minority-minority 
spin channels. With increasing bias voltage, the inelastic 
co-tunncling starts to enter the transport and leads to 
the decrease of x- Moreover, as the bias voltage reaches 
the threshold of sequential trans por t, the value of x a P~ 
proaches the lead-polarization p [26| , due to the appear- 
ance of sequential transport channels. 

The results of antiferromagnetic exchange interac- 
tion (Fig. 4) are different from the ferromagnetic case 
due to the two degenerate ground-states |1, 3/2; ±3/2) 
seen from Fig. 2(b), where the differential conduc- 
tance spectra have four evident peaks. With the in- 
crease of lead-polarization p (except p = 1) the po- 
sition of peak— 2 shifts toward the negative direction 
of the V-axis since the contribution is mainly from 
the transition of |0,2;±1) <^> |1, 5/2; ±3/2) instead of 
|0,2;±2) <s> |1, 5/2; ±5/2). The other three peaks 
(peak— 1, —3, —4) are resulted from the transitions of 
|0,2;±2) ^ |1, 3/2; ±3/2), |1, 5/2; ±5/2) & |2,2;±2), 
and |1, 3/2; ±3/2) |2, 2; ±2), respectively. For the case 
of p = 1, the position of peak-2 has a great change differ- 
ent from the peak-1 in the ferromagnetic coupling, where 
the state |0, 2; 2) does not participate in transition. Be- 
sides, I s and x exhibit similar behaviors with the ferro- 
magnetic case. 



B. Antiparallel configuration and current 
polarization reversal 

In the antiparallel configuration of magnetic leads pe- 
culiar behavior of the spin current is observed, which is 
shown in Fig. 5. From Fig. 5(a), it is seen that I c 
decreases monotonously from a maximum value to zero 
at a certain bias- voltage when p increases from to 1. 
Since transport occurs mainly through two spin channels 
(majority-minority and minority- majority) , for a large p 
the most probable transport process is that a spin-down 
electron of higher rate from the lead R tunnels through 
the molecule along with a spin-flip and molecular spin- 
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FIG. 5: (Color online) (a) Charge current I c , (b) differen- 
tial conductance G, (c) spin current I s , and (d) current po- 
larization x as a function of the bias voltage V for differ- 
ent lead-polarizations p with ferromagnetic exchange coupling 
(J = 0.2meV) and antiparallel lead-configuration. 



FIG. 6: ((Color online) (a) Charge current I c , (b) differen- 
tial conductance G, (c) spin current I 3 , and (d) current po- 
larization x as a function of the bias voltage V for different 
lead-polarizations p with antiferromagnetic exchange coupling 
(J = —0.2meV) and antiparallel lead-configuration. 



state lowering simultaneously. However, this process is 
limited by the lowest spin-state of molecule and thus the 
charge current finally vanishes [9j. Also, the G peaks in 
Fig. 5(b) are suppressed with the increase of p and the 
peak-2 appeared in Fig. 3(b) even disappears because the 
transition channels through states |0, 2; ±1) are blocked. 

The spin current I s (see Fig. 5(c)) shows a basin- 
like behavior in the region J-j- < I± similar to the QD 
case. The difference is that I 8 in the SMM can be sup- 
pressed by the lead-polarization p at some value, however 
it is enhanced always with the increase of p (p / 1) in 
the QD. Nevertheless, x displays a monotonic behavior 
with increase of p (p ^ 1) for a fixed bias voltage V 
seen from Fig. 5(d). It is found that the polarization 
X can have a high value at low lead-polarization p re- 
sulted from the exchange coupling between the electron- 
spin and the SMM. Moreover, in contrast to the parallel 
lead-configuration the current with \ — is contributed 
mainly from elastic co-tunneling processes around zero 
bias voltage, which belong to the majority-minority and 
minority-majority spin channels. With increasing V , the 
inelastic co-tunneling dominates the transport and x de- 
creases. Beyond the threshold value of sequential trans- 
port the peak-1 arises and x reaches the lowest value 
due to the transition |0,2;-2) |1, 5/2; -5/2). Then 
X increases since more transport channels enter the bias 
window. 

Fig. 6 is the I c — V diagram for antiferromagnetic 



exchange coupling, in which the I c curve (Fig. 6(a)) 
exhibits a similar behavior as in the ferromagnetic cou- 
pling case. The corresponding differential conductance 
G is shown in Fig. 6(b). Moreover, in the antiparal- 
lel lead-configuration the steady state probabilities with 
the negative eigenvalues of spin operator Sf ot become 
larger than that of positive eigenvalues resulted by the 
spin-flip process^, [TH, [TtJ • Therefore, the peak-1 and 
peak-3 of the G curve are mainly contributed by the 
spin-up electron transitions |0, 2;— 2) ^ 1,3/2;— 3/2) 
and |1, 5/2; -5/2) <S=> |2, 2; —2) respectively. On the 
other hand, the peak-2 (peak-4) is mainly induced by the 
transition |0, 2; —2) <S=> 1, 5/2; -5/2} (|1, 3/2; -3/2) <S=> 
2,2;— 2)) of the spin-down electron. The spin-current 
curve I s in Fig. 6(c) has a very interesting characteristic 
different from the ferromagnetic case that it possesses an 
alternate structure of two-plateau and one-basin. The 
first plateau appears between the positions of peak-1 
and peak-2, where the spin- up current is larger than the 
spin-down (positive x). In addition, the spin current 
shows a non-monotonic behavior with increasing lead- 
polarizations p, which is resulted from the competition 
between the spin-up and spin-down currents. Further- 
more, with increasing the bias voltage the basin emerges 
when the transition |0,2;— 2) |1, 5/2;— 5/2) starts to 
participate in the transport, in which the spin-down cur- 
rent becomes larger than the spin-up current (negative 
x). The second plateau of spin current appears as the 
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bias voltage increases to the value between the peak-3 
and peak-4. Finally when all transport channels open 
up, the spin current vanishes. In conclusion, the polar- 
ization of spin current x can be inverted by adjusting the 
bias voltage and the spin current exhibits rich multi-NDC 
behaviors }39| . The polarization reversal of spin current 
may have potential application in the spintronic device. 

IV. CONCLUSION 

In summary, the spin current through a SMM weakly 
coupled to two ferromagnetic leads is obtained in the 
nonlinear tunneling regime for both ferromagnetic and 
antiferromagnetic exchange-couplings. The complex en- 



ergy spectrum of the SMM results in rich properties of 
the currents and conductances compared with the sin- 
gle QD, which on the other hand provide magnetic sig- 
natures of the SMM. The most interesting observation 
is the spin-polarization reversal in the antiferromagnetic 
exchange-coupling case, which can be manipulated by the 
bias voltage. These theoretical results may be useful in 
the future development of spintronic devices. 
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